INTERNAL SHOCKS AND THE BLAZAR SEQUENCE 
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We consider the internal shock model as the dissipation mechanism responsible 
for the emission in blazars. It has been shown that this model is successful in 
reproducing the observed spectral energy distribution and the variability properties 
of a powerful blazar like 3C 279. However, the blazar family covers a wide range 
of spectral characteristics which appear to be correlated and the whole class can 
be seen as a sequence: the frequency and the intensity of the low energy versus 
high energy peak intensity increase with decreasing luminosity. We show that the 
internal shock model can satisfactorily account also for the properties of the low 
power blazars like BL Lac and Mkn 421 and it is successful in reproducing the 
blazar sequence. 



1 Introduction 

The discovery that blazars are strong 7-ray emitters together with the results 
from the multiwavelength campaigns of these sources have allowed to deepen 
our knowledge on these objects. We know that the Spectral Energy_Distribu- 
tion (SED) of blazars is characterized by two broad emission peaksEl strongly 
variable on different timescalesH. Two main radiation processes produce these 
peaks, i.e. the synchrotron at low frequencies and the inverse Compton at high 
energies (see e.g. □ for a review). The relative importance of the two peaks 
and theixJ-Ocation in frequency appear to be a function of the total power of 
blazars leading to a blazar sequence. The blazar emission is variable on 
different timescales, which are typically weekes-months for the radio emission 
of the order of a day for the 7-ray. 

Several studies, mainly based on the fitting of the observed SED, have al- 
lowed us to consistently derive the physical parameters in the emitting region 
of these sources, however many aspects in the understanding of relativistic 
jets remain open, most notably the jet energetics and particle acceleration. 
In order to explore these issues and their relationship we have considered a 



scenario where the plasma characteristics are not treated as free parameters, 
but are the result of the jet dynamics, thus relating the observed emission 
properties to the transport of energy along the jet. We have achieved that by 
quantitatively considering a dissipation mechanism within the jet responsible 
for the emission of blazars analogous to the standard mechanism proposed to 
explain gamma-ray burstsB: the internal shock model. This model was actu- 
ally originally proposed for blazars more than two decades agoB. The most 
important assumption of the model is that the central power engine produces 
energy which is channeled into jets in an intermittent way, though such a time 
dependent process cannot be easily inferred from first principles. The model 
assumes that different parts of the jet (shells), moving at different speeds can 
collide producing shocks, giving rise to the non-thermal radiation we see. This 
mechanism has a limited efficiency because only a small fraction of the bulk 
energy can be converted into radiation (unless tha-contrast between the two 
initial Lorents factor, is huge, see Beloborodov 20000; Guetta et al. 200ll3; but 
also Ghisellini 2002 □). In blazars, on the other hand, the radiative output has 
to be a small fraction of the energy transported by the jet (less than 10%), 
since most of the bulk energy has to go un-dissipated to the outer radio-lobes. 
The low efficiency required for blazars is a characteristic of the model, more- 
over the internal shock model can naturally explain other properties of these 
sources: a minimum dimension for the 7-ray source is required, as well as 
a minimum distance from the accretion disk, in order to avoid pair produc- 
tion. The internal shock model can naturally account for this, since the jet 
becomes radiative at (>10^^ — 10^^ cm) from the black hole. The observed 
large amplitude variability must be explained by a non steady state model: 
we believe that the shell-shell collision scenario is the simplest among them. 
The shells that have already collided once can collide again, at larger distances 
from the black hole with reduced efficiency, since the Lorentz factor contrast 
of the colliding shells decreases. This can explain why the luminosity of jets 
decrease with distance. Fig. 1 shows the efficiency as a function of distance 
for a simulation we have done for Mkn 421 E3. A detail study via numerical 

a' nulations of the predictions of this model in blazars has been carried out by 
for powerful objects (like 3C279). However the SED of this object does not 
represent the whole of the blazar family, which covers a wide range of spec- 
tral characteristics (luminosity, frequency of the peaks of emission and their 
relative intensity). We have applied our model to Mkn 421 and BL Lac itself 
which, besides being among the best studied BL Lac objects, are representative 
of the extremely blue and intermediate blazar class, respectively. Our aim is 
to see if the proposed scenario can explain the difference of the SED along the 
blazar sequence, more than to explain every detail of a particular spectrum of 
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a specific object. 

2 Internal shocks and spectrum 

2. 1 Generalities 

In this work we use the approximate model of the unsteady wind described 
in detail in Spada et al. (2001) and in the following we report only its most 
relevant assumptions. We consider a compact source, of typical dimension 
i?o ~ 10^^ cm, which produces an unstable relativistic wind characterized by 
an average luminosity Ljet- The emission from the jet is obtained by adding 
pulses radiated in a series of internal shocks that occur in the wind. The wind 
flow is described as a sequence of N = t^/ty shells, where is the overall 
duration of the emission of the wind and ty <C tw is the average interval 
between two consecutive shells. Each shell is characterized by four parameters 
(the subscript j denotes the jth shell): the ejection time tj, randomly extracted 
from a uniform distribution with an average value of ty ; the Lorentz factor Tj 
that is randomly extracted from a distribution between and Tm ; the mass 
Mj also randomly extracted; the width Aj ~ Rq. The dynamics of the wind 
expansion is characterized by a series of two-shell collisions in which the faster 
shells catch up with the slower ones in the outer parts of the ejecta. The 
energy released in each shock is assumed to be distributed among electrons 
and magnetic field with fractions eg and es, respectively. In this way we can 
determine the magnetic field strength needed to evaluate the synchrotron and 
the inverse Compton emission. Electrons are assumed to be accelerated at 
relativistic energies with a power law spectrum. 

2.2 Local spectra 

We refer the reader to Spada et al. (2001) for a detailed description of the 
derivation of the spectra that is obtained considering the synchrotron, the 
synchrotron self-Compton (SSC) and the external Compton radiation mecha- 
nisms. The only thing we would like to point out is that even in the case of low 
power BL Lacs we consider the presence of a source of soft photons external 
to the jet, which we identified with the emission reprocessed in the broad line 
region (BLR). We consider a luminosity L^xt = a-^disk (with a ^ 0.1) produced 
within -Rblr, corresponding to a radiation energy density (as measured inJie 
frame comoving with the shell) Ucxt = (17/12)aLdiskr^/(47ri?|Lj^c). (e.g. Ej). 
For simplicity, this seed photon component is considered to abruptly vanish 
beyond -Rblr- High and low power blazars have different line luminosities 
(and different photoionizing disk continua). According to Kaspi et al. (2000) 
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Figure 1; The result of a simulation of the internal shock model for Mkn 421. The average 
local radiative efficiency and the cumulative efficiency versus the distance from the black 
hole. The first "peak" in the efficiency corresponds to the first collisions between shells, 
while the "tail" at larger distances corresponds to second, third (and so on) collisions. The 
top of the plot shows a schematic representation of the jet, with grey levels proportional to 
the local efficiency. Prom Guetta et al. (2002), in prep. 
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, ^BLR and Ldisk are related by i?BLR « i^jgj^ with b ^ 0.7 This impUes 
that blazars with weaker broad emission hnes should have smaller broad line 
regions. In turn, this implies that, in low power BL Lacs, the first collisions 
between shells can occur preferentially outside i?BLR- 



3 Results and Discussion 

We applied our model to Mkn 421 and BL Lac itself. In Tab. 1 we list the input 
parameters used for the model, and for comparison we also report the same 
parameters for 3C 279. Via the numerical simulations the full time-resolved 
spectral behavior of the sources can be determined. The full time-dependent 
behavior of the simulated sources can be also examined through animations (see 



http://ares.merate. mi. astro. it/r^gabriele/421/index.htmi) in which the tempo- 



ral and spectral evolutions arc simultaneously shown. 

The results of this work show that the internal shock model can satisfac- 
torely account also for the properties of the low power blazars. We find that 
the key model parameters that need to be changed in order to reconstruct 
the sequence in the frame of the internal shock scenario, are indeed the broad 
lines intensity and the power carried in the jet which is proportional to the 
radiative luminosity. These parameters in turn regulate the SED shape, as 
they control the cooling efficiency of the emitting particles. The internal shock 
scenario in fact determines a characteristic time interval for the injection of 
relativistic electrons during each shell-shell collision, of the order of the dy- 
namical timescale t^ross ■ The intensity of the spectrum out of each collision is 
maximized at the end of this particular timescale. Two regimes are relevant: 
fast and slow cooling, corresponding to whether electrons of energy 7b {jbrneC^ 
is the minimum energy of the injected electrons) can (fast cooling) or not (slow 
cooling) radiatively cool in t cross- In highly powerful blazars the fast cooling 
regime applies in the inner regions (within the BLR) where also most of the 
power is dissipated. Consequently the peak frequencies are produced by the 
electrons of energy 7t,meC^ (as even these electrons can cool in tcross)- In the 
lowest powerful blazars instead the slow cooling regime applies to all collisions. 
Consequently only the electrons with the highest energy can cool in tcross (even 
for the most powerful collisions): the peak frequencies thus shift to very high 
values. Between these two extremes there is the possibility of BL Lac objects 
with broad lines of intermediate intensity originating at a distance from the 
accretion disc within which very few (but not zero) shell-shell collisions can 
occasionally take place. This is the case of BL Lac itself and all the blazars of 
intermediate luminosity. Observationally this corresponds to a SED produced 
by internal shocks of moderate Compton to synchrotron luminosity ratio as 
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Table 1: Input parameters of the model for the 3 blazars. 



the seed photons for Compton scattering are provided only by the synchrotron 
photons (colhsions outside the BLR). The rare coUisions within the BLR will 
give rise to a dramatic change in the SED characterized by a large increase of 
the Compton component as in the case of BL Lac itself. 

One can conclude that the internal shock scenario explains the blazar sequence 
and the main characteristics of the blazar behaviour. While observationally we 
were able to see that the important quantity is the ratio between the observed 
luminosity in the jet and the disc (thus BLR), the internal shock scenario al- 
lows to i) directly connect the radiated jet luminosity with the effective power 
carried by the jet and ii) account for the preferred distance at which most of 
the jet luminosity is produced. While the latter distance is qualitatively the 
same for high and low power blazars, the BLR is instead located at very differ- 
ent distances in the two sub-classes of objects, as determined by the ionizing 
luminosity. 
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